Formation of fermionic molecules via interisotope Feshbach resonances 
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We perform an analysis of recent experimental measurements and improve the lithium interaction 
potentials. For ®Li a consistent description can be given. We discuss theoretical uncertainties for the 
position of the wide ^Li Feshbach resonance, and we present an analytic scattering model for this 
resonance, based on the inclusion of a field-dependent virtual open-channel state. We predict new 
Feshbach resonances for the ®Li-^Li system, and their importance for different types of crossover 
superfluidity models is discussed. 
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Resonances in cold atomic gases offer the key to con- 
nections with challenging condensed matter physics. In 
particular, resonances make lithium atomic systems very 
versatile. The first Bose-Einstein condensates (BEG) 
were quite small in number due to a negative scatter- 
ing length [l|, later Feshbach resonances Q] have been 
used to create condensates with positive scattering length 
and to generate bright solitons 0j 01 by changing the 
scattering length a back to negative. Even more in- 
teresting is the usability to form molecules, since the 
Feshbach resonance results from bringing a molecular 
state on threshold. The connection between fermionic 
atoms and composite bosons (molecules) has great im- 
pact for the study of the well-known crossover prob- 
lem between BEG and Bardeen-Goopcr-Schrieffer (BGS)- 
type superfluidity[l,|g,[3. 

In this paper, we study several Feshbach resonances 
in the lithium system. We first review the knowledge 
of the interatomic interaction potentials, and use exper- 
imental data as input to improve these potentials. We 
discuss the special situation for the wide ^Li resonance, 
where the background scattering length depends strongly 
on the magnetic field. This will be interpreted as a field- 
dependent virtual state (a second resonance) which is 
situated close to threshold. The full energy-dependent 
scattering process can be parametrized according to a 
simple analytical model that encapsulates both field- 
dependent resonances. Further, we apply our knowledge 
of the lithium interactions to a system of ^Li-^Li, and 
find several Feshbach resonances which are accessible in 
current experimental setups. The underlying molecular 
state is of a composite fermionic nature, which allows for 
a new type of crossover physics — the transition of an 
atomic BEG to a molecular Fermi-type of superfiuidity 
Recently, Feshbach resonances in a heteronuclear Bose- 
Fermi mixture have been observed [a, |9j, where polar 
fermionic molecules underly the resonance state. 

For an accurate prediction of resonance properties, 
we need a detailed understanding of the actual inter- 
atomic potentials. Here we describe how we improved 
the precision of existing potentials by using recent ex- 
perimental measurements as input. The potentials can 
be divided in two radial intervals. For large interatomic 
separations r the potential is given by the sum of the 
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FIG. 1: (color online) Dependency of the wide-resonance po- 
sition Bo on Ce, according to the second analysis (see text). 
The shaded area indicates the associated uncertainty. The 
curve indicated with A (B) is obtained making use of the 
exchange energy according to [ill (O)- 



dispersive Van der Waals tail V^'^^{r) 
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Cs/r° —Cm/r and the exchange contribution Vg^{r) = 
(-l)^+iC"=^r^/2"-ie-2«'- [l3, resulting in two poten- 
tials: a singlet (5" = 0) and a triplet {S = 1) p otential. 
The coefficient C"^^ is taken from Refs. [13,^3, a is di- 
rectly related to the ionization energy a^/2 [l^l, and Cg 
and Cio are taken from Refs. [lj|. For smaller r we use 
the model singlet and triplet potentials which have also 
been used in Refs. [l5l[la|. 

These short-range and long-range potentials are 
smoothly connected at r = ISag, with oq the Bohr radius. 
To overcome the inaccuracies of the short range poten- 
tials, we make use of the accumulated phase method [l5| . 
A boundary condition is applied on the partial-wave ra- 
dial wave functions at r = 17.5ao in the form of a WKB 
phase (/)5,t(^,^) = 0s^t(S,£) + A0s,t- The first term on 
the right is calculated by radial integration of the model 
potential up to 17.5ao and is expected to account for the 
energy and angular momentum dependence of the accu- 
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FIG. 2: (color online) Energy of the dressed (quasi-)niolecular 
state as obtained with a full coupled-channels calculation 
(dots), with respect to the |/i,m/j) g) |/2,m./2) = 1/2, 1/2) 8 
1/2, —1/2) threshold. The straight line is the bare closed- 
channel energy eb{B), and the energy of the dressed state 
according to our model is indicated by the curved line. Inset: 
corresponding background value of the scattering length as a 
function of the magnetic field B. 



mulated phase to a sufficient degree of accuracy. The 
second term is an energy and angular momentum inde- 
pendent shift of the phase, determined from experimen- 
tal data. These corrections A(J>s,t to the accumulated 
singlet and triplet phases can be converted to the more 
physical quantities vds.dt, which arc the fractional vi- 
brational quantum numbers at dissociation. 

We determine the free parameters of our interaction 
potentials vds^ ^dt, and Cg from experimental input 
by means of a y^ minimization. An inter-isotope anal- 
ysis in which ''Li is related to ^Li by means of a simple 
mass-scaling relation failed, yielding inconsistent results 
for vds- This is a strong indication of a breakdown of 
the Born-Oppenheimer approximation for the singlet po- 
tential. Such a breakdown was demonstrated in detailed 
spectroscopy |l3 . We therefore avoid mass scaling of the 
singlet potential, and we perform two different analyses. 
In the first analysis, we only take ^Li data into account. 
In the second analysis we investigate ^Li as well, how- 
ever, we only do a mass scaling for the triplet potential. 
Our total set of ^Li experiments comprises 6 data points. 
The zero-crossing of the scatter ing length of a system in 
the two lowest hyperfine states [l^, [l3 ■ In the same spin 
state configuration, the positions of the narrow [23 and 
wide Q Feshbach resonances. The measurement of the 
scattering length in the lowest and third to lowest hyper- 
fine state |2ll | and the binding energy of the most weakly 
bound triplet state pa |. 

In our first analysis we obtain a minimum in the re- 
duced x^ distribution of x^ = 0-5. The corresponding 



parameter values are vj^s = 0.3496(5), equivalent with 
a singlet scattering length as ~ 45.3(l)ao and vdt = 
0.9954(2), corresponding to ar = -2025(70)ao. For 
the leading dispersion coefficient we find Cg = 1388(6) 
at. units. This result has been obtained with a C^^ of 
Rcf. O. When we weaken C"=^ to the value of Rcf. [ll|, 
we find that the optimal Cg is shifted to Cg = 1390 
at. units. The scattering lengths found are consistent 
with previous determinations [l^, |23|- Our Cg coeffi- 
cient agrees with the values found in ah initio calcula- 
tions [mIII. 

The objective of a second analysis [2J| is to reevalu- 
ate the position of the wide s-wave Feshbach resonance 
of ^Li in the two lowest hyperfine states, without mak- 
ing use of the experimental result i3o=822G |3|. Hereto 
we want to combine all available cold collision data on 
lithium, and we add the positions of three p-wave Fesh- 
bach resonances, which have been measured recently |25j . 
and experimental data of ^Li [3, [23 , to the set of exper- 
iments. In this combined isotope analysis we perform a 
mass scaling procedure for the triplet boundary condi- 
tion only. As explained above, we will not mass-scale the 
singlet potential but rather optimize the boundary con- 
ditions for the singlet potential independently for the two 
isotopes, by making vds,7 a free parameter independent 

Since the p-wave resonances are measured with high 
accuracy, we also allow for small corrections of the angu- 
lar momentum dependence of (j)^{E,£) via the parame- 
ter A(/)g by means of an addition Acj)g ■ l{l -f 1) cf. |27l| . 
By optimizing the interaction parameters {vosfi^ ^ds,7i 
V]jT.&; A0s) ^'^^ various fixed values of Cg we are able 
to obtain a minimum reduced x^ of 0.7. The depen- 
dency of x^ on Cq is rather weak and therefore the set 
of experiments does not restrict the Cg coefficient to an 
acceptable degree. However, the minimal x^ occurs for 
Cg = 1390.6 at. units close to ah initio values, when using 
(jex f].Qj-Q ^Qf (n|, positioning the wide Feshbach reso- 
nance at Bq = 826 G. We estimate 808 G < Bq < 846 G 
for 1388 < Cg < 1393 at. units, see Fig. HI For C^^ from 
Ref. [13I we find a minimum x^ = 0.5 for Cg = 1395.6. 

From now on, the properties of the Feshbach reso- 
nances will be derived from our first analysis. First, we 
study the wide Bq = 822 G Feshbach resonance in the 
two lowest hyperfine states of ^Li. This resonance is quite 
remarkable for two reasons: it has a large width of the 
order of 100 G, and its background scattering length Obg 
is strongly depending on the magnetic field, which can be 
seen from the inset of Fig. [3 At zero field, Obg is 3 and 
positive, while for large field values Obg is large and neg- 
ative, indicating the presence of a nearby virtual state in 
the open-channel subspace P |23|. Consequently several 
important quantities, such as the S and T matrices which 
summarize the collision process, depend non-trivially on 
the collision energy i?. We will here apply the model 
discussed in Ref. [23| to this s-wave Feshbach resonance. 
This model takes the virtual state into account explic- 
itly, and gives an analytical description of all important 



two-body quantities near the Feshbach resonance. 

In general, the relation between the background scat- 
tering length Obg, the range of the potential a^^, and the 
virtual state pole k^s, is given by Obg = a^„ — l/^vs- The 
range of the potential is related to the Van der Waals co- 
efficient Cg, and does not depend on the magnetic field. 
Therefore, we account for the field dependence of Obg 
by generalizing the model of Ref. [231 to the case of a 
field-dependent virtual-state Kvs(-S)- The complex en- 
ergy shift is then given by 
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where Avs{B) is related to the coupling matrix element 
between the open-channel virtual state and the closed- 
channel bound state responsible for the Feshbach reso- 
nance. Our wavenumber units are such that E = k'^. 
The total scattering length is then given by 
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where eb{B) = Afj.'^'^s{B — Bq) is the energy of the bare 
closed-channel bound state, A/i™''^ = 2.0^b the mag- 
netic moment, hb the Bohr magneton, and Bq ~ 539.5 G 
the field where the bare closed-channel energy crosses 
threshold. For a fixed B value, the parameters ^vs(^) 
and Kys{B) are obtained by fitting Eq. Q to the coupled- 
channels result for a{B), using two close-lying data points 
where we assume A^g and Kvs locally constant. Repeat- 
ing this for every field value, we obtain explicit expres- 
sions for the non-trivial energy dependence of the com- 
plex energy shift, using only the zero-energy information 
contained in a{B). Note that there are only two free 
parameters, and ahg{B) is fixed once Kvs(-B) is known. 

The fit functions are summarized in Table Q] In Fig. 13 
we compare the dressed (quasi-)molecular state calcu- 
lated by coupled-channels methods and by our analyti- 
cal model, which agrees excellent. Therefore this model 
can be used to analytically describe the S and T ma- 
trices (2^, scattering phase-shifts, etc., with similar pre- 
cision as a full coupled channels calculation, for a large 
range of energies and magnetic fields. We note that also 
the narrow 543 G Feshbach resonance (width of order 
0.1 G) can be described by Feshbach theory. Here, how- 
ever, an easier description is possible based on only one 
single background part and a single resonance state, since 
the narrow resonance has a 'local' background scattering 
length of the order of a^ . For narrow resonances |29ll30l |. 
the typical resonance features in the continuum scatter- 
ing which depend on the details of the potential, are very 
important to the BEC-BCS crossover physics, since they 
are visible for energies less than the Fermi energy, result- 
ing in a non-universal crossover picture. 

Now we turn to the final topic of this paper, the study 
of Feshbach resonances in a ^Li fermion - ^Li boson mix- 
ture. We use again our first analysis of the interaction 
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FIG. 3: (color online) Two thresholds for the ^Li-^Li system 
which show the coincidence of two Feshbach resonances. In- 
set: scattering length as a function of magnetic field, for the 
two coinciding Feshbach resonances. 



TABLE I: Parameters of the virtual-state model. For the 
magnetic fields of interest (530 G< B < 830 G), the field 
dependent parameters are given by a third-order polynomial 
fit, co + ciB + C2B^ + C3B^. For Avb(B) the units are given by 
[c„] — K^G"", and for k.vb{B) the units are given by [c„] — 
K^'^G~". The range of the potential is given by a^^ = 45 ao. 
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parameters for ^Li, and perform a simple mass scaling of 
the accumulated phases for the ^Li-'^Li system. As dis- 
cussed before, this leads to less accurate predictions for 
the resonance positions as suggested by the accuracies 
from the ^Li interaction parameters. We estimate the 
inaccuracies in the mixed isotope resonance positions to 
be of order a few Gauss, due to inaccuracies of the mass- 
scaling relations. 

Feshbach resonances in the ^Li-^Li system have been 
studied before in Ref. [ig , where only the case of magnet- 
ically trappable atoms was investigated. Moreover, those 
resonances are accompanied by large inelastic exchange 
losses. We investigate only ^Li-^Li hyperfine state com- 
binations where exchange losses are absent. Within these 
boundaries, there are still numerous resonances present, 
and we restrict ourselves to the most interesting results. 
The *^Li-'^Li |l/2, 1/2) |1, 1) channel has the lowest en- 
ergy in the two-body hyperfine diagram. Therefore, this 
channel will not suffer from magnetic dipolar relaxation. 
We find five Feshbach resonances at magnetic field val- 
ues of 218 G, 230 G, 251 G, 551 G, and 559 G. Mea- 
surements of these resonances might provide the missing 



information to exactly locate the position of the wide ^Li 
resonance, as all mixed resonances arise from the same 
underlying bound state in the triplet potential. 

Feshbach resonances between bosons and fermions give 
rise to a fundamentally different type of crossover physics. 
Already some work on the interactions in Bose-Fermi 
mixtures can be found in the literature pUll32l|. however, 
in order to describe correctly the many-body physics and 
interactions close to resonance, more research is needed. 
To make the system feasible, some requirements have 
to be fulfilled. For a stable BEC, a positive scattering 
length for the bosons is required. Approaching the res- 
onance from the atomic side, the mixed boson-fermion 
scattering length will become negative, and stability of 
the system could become and issue. At the other side 
of the resonance, it is important that the effective in- 
teraction between the fermionic molecules is attractive. 
Point-like composite fermions do not undergo s-wave col- 
lisions. However, close to resonance the molecules are 
long stretched, and an effective interaction mediated via 
the bosons could be possible. An approach similar to 
Ref. [33 could be conclusive on this and on the expected 
dependence of the inelastic rate coefficient on the scatter- 
ing length. Also, the effect of Pauli blocking will not be 
as strong as in Ref. [23 since a three-body decay process 
with only one fermion and two bosons involved is possi- 
ble, however, still a reduction with respect to the pure 
bosonic case could be expected. 

Another interesting situation occurs when two Fesh- 
bach resonances coincide. Further away from resonance, 
where the size of the molecule is comparable to the size 
of the potential, s-wave collisions arc not allowed for 



these composite fermions. Therefore, in order to pre- 
serve superfluid behavior in this region, two different 
molecular spin configurations are needed to allow for s- 
wave collisions between molecules. Moreover, it is desir- 
able that the two resonances responsible for the molecule 
formation coincide. Such a coincidence can be found 
from Fig. O where the ^Li-'^Li |l/2, -1/2} ® |1,1} and 
3/2, —3/2) (g) |1, 1} bound states and thresholds are plot- 
ted as a function of magnetic field. Every crossing of a 
bound state with threshold indicates the position of a 
Feshbach resonance. It can be seen that a.t B = 305 G 
two Feshbach resonances coincide. These two threshold 
channels have the same bosonic ''Li state, but a different 
fermionic ^Li state. This coincidence is systematic, and 
will not depend on the interaction parameters. The scat- 
tering lengths as a function of the magnetic field can be 
seen in the inset of Fig. |3| 

In conclusion, we analysed recent experimental mea- 
surements for lithium, and showed that mass scaling be- 
tween ^Li and ^Li fails for the singlet potential. We inves- 
tigated uncertainties in the wide ^Li resonance position, 
and demonstated an analytical model for this resonance 
that includes the nearby virtual state. Finally, wc showed 
that ^Li-^Li mixtures feature accessible Feshbach reso- 
nances, giving rise to fermionic molecules, yielding new 
BCS-BEC crossover physics. 
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